Limits on Low Energy Photon-Photon Scattering from an 
Experiment on Magnetic Vacuum Birefringence 



M. Bregant,^ G. Cantatore/ S. Carusotto,^ R. Cimino,^ F. Delia Valle/ 
G. Di Domenico,^ U. Gastaldi,^ M. Kamza,^ V. Lozza/ E. Milotti/ 
E. Polacco,^ G. Raiteri/ G. Ruoso,^ E. Zavattini*/ and G. Zavattini^ 

g (PVLAS Collaboration) 



O 
<N 



^INFN, sezione di Trieste and Dipartimento di Fisica, 



^ Universitd di Trieste, Via Valeria 2, 34127 Trieste, Italy 

^ ^INFN, Sezione di Pisa and Dipartimento di Fisica, 

O Universita di Pisa, Via Buonarroti 2, 56100 Pisa Italy 

<N 

^INFN, Laboratori Nazionali di Frascati, 

Via E. Fermi 40, 00044 Frascati, Italy 

Ph ^INFN, Sezione di Ferrara and Dipartimento di Fisica, 

r^. Universitd di Ferrara, Polo Scientifico, 

^ Via Saragat 1 C, 44100 Ferrara, Italy 

> 

'^INFN, Laboratori Nazionali di Legnaro, 

cn 

O viale dell 'Universitd 2, 35020 Legnaro 

cn 

lyr^ ^INFN, sezione di Trieste, Via Valerio 2, 34127 Trieste, Italy 

O 

OO (Dated: May 20, 2008) 

o 



Deceased January 9, 2007 



1 



Abstract 

Experimental bounds on induced vacuum magnetic birefringence can be used to improve present 
photon-photon scattering Umits in the electronvolt energy range. Measurements with the PVLAS 
apparatus (E. Zavattini et al, Phys. Rev. D 77 (2008) 032006) at both A = 1064 nm and 
532 nm lead to bounds on the parameter Ag, describing non linear effects in QED, of yli^^^^-* < 
6.6 • 10-21 @ 1064 nm and ^i^^^^ < 6.3 • lO^^^ T'"^ @ 532 nm, respectively, at 95% confidence 
level, compared to the predicted value of Ag = 1.32-10-2^ T^^. The total photon-photon scattering 
cross section may also be expressed in terms of A^, setting bounds for unpolarized light of a^^''^^^ < 
4.6 • 10"^^ m^ and a^^^^^ < 2.7 • 10~^^ m^. Compared to the expected QED scattering cross section 
these results are a factor of ~ 2 • 10^ higher and represent an improvement of a factor about 500 on 
previous bounds based on ellipticity measurements and of a factor of about 10^*^ on bounds based 
on direct stimulated scattering measurements. 

PACS numbers: 12.20.Fv, 07.60.Fs, 14.80.Mz 
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I. INTRODUCTION 



Classical electrodynamics in vacuum is a linear theory and does not foresee photon-photon 
scattering or other non linear effects between electromagnetic fields. Before quantum electro- 
dynamics (QED) was formally complete, Euler and Heisenberg, and Weisskopf, realized that 
vacuum fluctuations, permitted by the uncertainty principle, lead to non linear effects: 4 pho- 
tons can couple via fermion loops. Such non linear effects were first calculated in 1936 [Ij and 
can be described by an effective Lagrangian Lehw which, for field strengths well below their 
critical values {B < ficrit = m^c'^/eh = 4.4 ■ 10^ T, E < E^nt = m^c^/eh = 1.3 • 10^^ V/m) 
and for photon energies below the electron mass, can be written as (in S.I. units): 

2 



^EHW — 




(1) 



= = 1.32 ■ 10-24 (2) 



where the parameter describing the non linearity is 

2 a^Xl 

with Ae being the Compton wavelength of the electron, a = / {hcATreo) the fine struc- 
ture constant, irif. the electron mass, c the speed of light in vacuum and /io the magnetic 
permeability of vacuum. 

Maxwell's equations are still valid provided the constitutive equations are applied to the 
total Lagrangian density L = Lciass + -Z^ehw to derive the displacement vector D and the 
magnetic field intensity vector H 

Co dE 

H . (3) 

One of the yet to be measured effects predicted by the Lehw correction is that vacuum 
will become birefringent in the presence of an external magnetic and/or electric field. For 
example, in the case of a beam propagating perpendicularly to an external magnetic field, 
if ny and n± indicate the index of refraction for polarizations respectively parallel and 
perpendicular to the field direction, the birefringence can be expressed as [21 13] 



nil 



= An(Q^D) = 3A,B^ (4) 
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which is extremely smaU: with a field intensity of 5 T, An*^^™) ^ 10~^^. Another process 
described by the same Feynman diagrams as magnetically induced vacuum birefringence is 
photon-photon scattering. Figure [T] shows the Feynman diagrams for both photon-photon 
scattering and field induced vacuum birefringence. In general, the connection between the 



Figure 1: First order Feynman diagrams for both photon-photon scattering and magnetically (elec- 
trically) induced vacuum birefringence. 

index of refraction n of a medium and the photon scattering amplitude in the forward 
direction for photons with energy E^, /(-i? = 0, E^), is (see for example 

n = l + ^NfiO,E,) (5) 

where N is the average number density of centers of scattering and k is the photon wave 
number. Applied to photon-photon scattering of linearly polarized photons, the center of 
mass forward scattering amplitude of ingoing and outgoing photons all having parallel polar- 
izations, /|j^^°'*(0, E'^), and the one in which the two incoming photons have perpendicular 
polarizations as do the ongoing photons, f^^^\o, E^) are, respectively [S] 

p{QED)/^ . _ 32 a% / -^7 V _ 16/io A 3 



/p^""^(0,E,) = [—2]= i^2^eE,' (6) 

KQED).^ X _ 56 a% / ^7 _ . 3 



where it is apparent that the scattering amplitude is proportional to A^,. The authors of [5] 
also show that is proportional to the energy density of the scatterer field (electric and/or 
magnetic) and inversely proportional to the photon energy in the center of mass reference 
frame. From the scattering amplitude one can find the differential cross section 

^(^,i?,) = |/(^,i?,)P (8) 
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and the total cross section which depends on A^. For unpolarized hght one finds [HI [TJ [U [9l [TO] 



77 v"7/ 452 57r" V^ecV 207r h^d' 

The connection between the total photon-photon cross section and vacuum birefringence, 
hence the parameter describing non linear QED effects, makes non linear QED searches 
via ellipsometric techniques very attractive. Limits on A^ from ellipsometric data can there- 
fore be directly translated into photon-photon scattering limits. 

It is interesting to note that in a post-Maxwellian framework [H] the Lagrangian den- 
sity LpM describing nonlinear electrodynamic effects in vacuum is parameterized by three 
parameters ^, rji and 772: 

In this parameterization ^ = 1 / B^^.^^ , and rji and t]2 are dimensionless parameters depending 
on the chosen model. In the Euler-Heisenberg electrodynamics rj^^^^ = ^rj^^^^ = a/(457r), 
a being the fine structure constant. 

By substituting the post-Maxwellian generalization into equations ([s]) one finds that the 
birefringence induced by a transverse magnetic field is (to be compared with equation ^ 



An(P^) = 2e(r/2 - r7i)52 (11) 
whereas the forward scattering amplitudes given in expressions (6) and (7) become 

fr\o,E,) = ^,(mE,:' (13) 

Birefringence is therefore only sensitive to the difference 772 — Vi whereas the two forward 
scattering amplitudes fl-^^\o, E^) and f^^\o, E^) are proportional respectively to rji and 
?72- At scattering angles different from = it remains true that fl^^\'d, E^) is propor- 
tional to 7]i but f^^\i!}, Ey) will now depend on a combination of rji and 772 which never 
cancels. Therefore, for example, in the Born-Infeld model [12] where rji = ri2, magnetically 
induced birefringence is not expected even though photon-photon scattering is. Although 
very promising for detecting nonlinear electrodynamic effects, the ellipsometric technique 
alone it is not sufficient to determine the two independent quantities ^rji and ^772- On the 



other hand direct photon-photon scattering with defined polarization states can. It is clear 
how both techniques are complementary. 

Assuming the Euler-Heisenberg Lagrangian density, in this paper we will present the best 
limits on a^^ at low energy available today. 



II. APPARATUS AND METHOD 



The general scheme of a sensitive ellipsometer searching for magnetically induced bire- 
fringence is presented in Figure [2j A polarizer defines the polarization of the beam, of power 
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Figure 2: Schematic layout of a sensitive ellipsometer. See text for description. 



Jin, before it enters the magnetic field region where it acquires an ellipticity ip. The ellipticity 
is made time dependent by modulating the magnetic field with angular frequency ^Mag (see 
text below). Two mirrors compose either a multi-pass or a Fabry-Perot cavity to increase 
the optical path within the magnetic field region. The beam then passes first through a 
modulator, where it acquires a known ellipticity ^ modulated at frequency c^Mod, and then 
through an analyzer. The transmitted power /out is then detected and analysed. 



A. Heterodyne technique 

For the purpose of our discussion let a laser beam propagate along the Z axis and let the 
incoming (linear) polarization define the X axis (Figure [s]). Considering the coherence of 
our light source, the Jones matrix formalism will be used. The Jones matrix for a uniaxial 
birefringent element is given by 

(1+2-0 COS 2^ lib sin 2'd 
itjj sin 2'd 1 — i-ip cos 2i!} 




6 



X 



X 



n. 



'II 




Birefringent 



y 



Z 



Y 



Figure 3: Reference frame for the calculations using the Jones matrix formalism. 



where represents the angle between the slow axis {n\\ > n±) of the medium and the X 
axis and <S 1 is the induced ellipticity acquired by the light, given by 



^^n^^^MTH::^ (15) 

2 A 



with — 9?^ the phase delay between the parallel and perpendicular polarization components 
acquired in a length L. 

Given an input beam whose electric field after the entrance polarizer is = -E'in(o) the 
electric field after the birefringent medium will be 



Assuming no losses, the power /out after the analyzer (polarizer crossed with respect to the 
entrance polarizer) will therefore be 



The power is proportional to iJj'^ and, whether i) is constant in time or not, results in an 
unmeasurably small intensity component. 

To linearize the term proportional to the ellipticity signal ip to be detected, one can add 
a known time varying ellipticity (^{t) using an ellipticity modulator. The Jones matrix for 
the modulator is the same as BF set at an angle of7r/4('0<^C<Cl): 




I out = -^in \itpsm2'df 



(16) 




(17) 
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and the resulting vector describing the electric field after the modulator will be 

^0 = ^in-MOD • BF • Q (18) 

f l + ^V'COs21?-V'C(^)sin2^? \ ^^^^ 
I <(t) + lip sin 2i& - C{t)ip cos 2i9 j 

Neglecting second order terms, the power /out after the analyzer will be 

/out(t) = /in |<(t) + #sin2i?|' ~ Ji, [C{tf + 2C(t)V^sin2i?] (20) 

which now depends linearly on the ellipticity ■?/'■ To complete the discussion, one finds ex- 
perimentally that static and slowly varying cllipticitics, in the following indicated as oi{t), 
arc always present in an actual apparatus and that two crossed polarizers have an intrinsic 
extinction ratio cx^, mainly due to imperfections in the crystals they are made of. Further- 
more, losses in the system reduce the total light reaching the analyzer. Therefore, taking 
into account an additional spurious ellipticity term a{t) (since a,^, C ^ 1 these terms com- 
mute and therefore add up algebraically) and a term proportional to cr^, the total power at 
the output of the analyzer will be 

Iont{t) = h [(7^ + \iC{t) + ^V' sin 2'& + ia{t)\^] ~ 

~ lo [a^ + C{t? + a{tf + 2C{t)iP sin 2^ + 2C{t)a{t)] (21) 

where Iq represents the power of light reaching the analyser. 

To be able to distinguish the term ((t)a(t), which is usually largest at low frequencies, 
from the term ({t)ilj sin 2'&, the term of interest ip sin 2'& is also made to vary in time. This can 
be done by either ramping the magnetic field intensity (varying therefore ip) or by rotating 
the magnetic field direction (varying -d). The final expression, explicitly indicating the time 
dependence of ip and for the power at the output of the analyzer is therefore 

/out(t) = lo [a'' + C{tf + a{tf + 2C(t)V^(t) sin 2^{t) + 2C,{t)a{t)\ (22) 
B. Optical path multiplier 

To further increase the ellipticity induced by the birefringent region one can increase the 
number of passes through it. Either a multi-pass cavity or a Fabry-Perot cavity can be used 
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for this purpose. In the PVLAS experiment described below, a Fabry-Perot has been chosen. 
In a multi-pass cavity the induced ellipticity is proportional to the number of passes A^pass 
through the region. With a Fabry-Perot cavity the calculation is not immediate since one 
is dealing with a standing wave. 

Let t, r be the transmission and reflection coefficients, and p the losses of the mirrors of 
the cavity such that + + p = 1. Let d be the length of the cavity and 6 = Aird/X the 
roundtrip phase for a beam of wavelength A. Then the Jones matrix for the elements of the 
ellipsometer after the entrance polarizer is 

oo 

ELL = A ■ SP ■ MOD ■ t'^e'^^^ ^ [BFV^e*''] " • BF (23) 

n=0 

(oo\ 

where A = is the analyzer Jones matrix and SP describes the spurious ellipticity. 

Because < 1, ELL can be rewritten as 

ELL = A ■ SP ■ MOD ■ t^e'^^^ [l - BFV^e*^] ■ BF (24) 

with I the identity matrix. With the laser phase locked to the cavity so that 6 = 27rm, 
where m is an integer number, the electric field at the output of the system will be 



E,^, = Ei,-ELL ■ ( = ^inT^ f ° 2 I (25) 

/ t^+P \ ta{t) + <(t) + sin 2^ 



and the power, including losses, 

2 



ia{t) + <(t) + I [ \ ^ ^ ] i> sin 2d 



(26) 



This expression is at the basis of the ellipsometer in the PVLAS apparatus. Small elliptici- 
ties add up algebraically and the Fabry-Perot multiplies the single pass ellipticity ip sin 2d, 
generated within the cavity, by a factor (1 + r^)/(l — r^) 2jF/7r, where is the finesse 
of the cavity. The ellipticity signal to be detected is therefore ^ = (2jF/7r)?/' sin2'(9. Typical 
values for the finesse of the PVLAS cavity are ~ 10^. 

In the PVLAS experiment, ((t) = Co cos(co'Mod^ + Ouod) and the magnetic field direction 
is rotated at an angular velocity ^uag- A Fourier analysis of the power /out(^) of equation 



(26) results in four main frequency components each with a definite amplitude and phase. 



These are reported in table |T| 



Table I: Intensity of the frequency components of the signal after the analyzer A. 



Frequency 



Fourier component 



Intensity //q 



Phase 



DC 

^Mod 
2a;Mod 



I 



<^Mod 



'^'^Modi2r2Ma 



2aDcCo 

C0V2 



t'Mod 
^'Mod ± 2^Ma 
2^Mod 



The presence of a component at c<jMod±2f2Mag in the signal identifies an induced elUpticity 
within the Fabry-Perot cavity. Furthermore the phase of this component must satisfy the 
value in table HI 



C. PVLAS Apparatus 
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Figure 4: Schematic layout of the PVLAS apparatus. See text for description. 

A description of the PVLAS apparatus, shown schematically in Figure |4], can be found 
in [131 HH US] • The magnetic field is provided by a superconducting dipole magnet which is 
placed vertically and rotates around its axis, at a typical frequency of 0.3 Hz. The magnetic 
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field therefore lies in the horizontal plane, the field region is 1 m long, and the maximum 
field intensity is 5.5 T. The ellipsometer develops vertically: polarizer and entrance cavity 
mirror are supported by the lower optical bench, whereas output cavity mirror, modulator 
and analyzer are in a vacuum chamber on the upper optical bench. The lower optical 
bench is in a pit whose floor is a concrete slab resting on four 14 m long pillars buried 
in the ground. The slab and pillars are therefore seismically isolated with respect to the 
surrounding hall floor and building. The upper optical bench is sustained by a granite tower 
7 m high also standing on the concrete slab. The upper and lower vacuum chambers are 
connected by a quartz tube 2.5 cm in diameter which passes through the warm bore of the 
cryostat containing the magnet. The magnet and turntable are supported by a concrete 
beam crossing over the pit and resting on the hall floor. Thus, mechanical vibrations due to 
the rotating magnet reaching the optical system will be greatly suppressed and should not 
cause excess elhpticity noise. 

The vacuum system is based on two liquid N2 traps combined with Ti sublimation getter, 
and pressure is kept at the level of P ~ 10"^ mbar during measurements. For test purposes 
the vacuum chamber can be filled with gases at known pressure, measured with a set of 
capacitive transducers. The presence of the gas gives rise to a known magnetic birefringence 
via the Cotton Mouton effect [T6l[T71[T8]. 

The laser source is frequency locked to the Fabry-Perot cavity using a modified Pound- 
Drever-Hall technique [12]. Two different light sources were alternatively used: an infrared 
Nd:YAG laser emitting 800 mW at 1064 nm (infrared), and its frequency doubled secondary 
output of 80 mW at 532 nm (green). The cavity parameters were as follows: finesse J-'ioga = 
70000, output power Pio64= 60 mW for the infrared and JF532 = 37000, P532=1.5 mW for 
the green. 

The light transmitted by the analyzer is detected by a photodiode connected to a low 
noise current amplifier and the signal is then sent to both a spectrum analyzer, for online 
monitoring of the apparatus, and to a lock-in amplifier demodulated at c^Mod- To make 
the analysis independent from the instability of the rotation frequency of the turntable 
sustaining the magnet, the table perimeter is equipped with 32 equally spaced trigger marks. 
The output of the lock-in amplifier is acquired at the passage of each trigger mark, therefore 
maintaining the coherence of the searched signal even for long integration times. 
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D. Noise considerations 



In the presence of a signal above background with the correct Fourier phase, the elhpticity 
\1/ = (2jF/7r)?/' can be calculated from Iq, from the Fourier components /ajMod±20Mag' from 
hujMod the average of the two sideband signals: 

^ 1 I -^'^Mod+2f!Mag _|_ -^I^Mod — 25lMag j (27) 

2 \ A/2/o/2c^Mod A/2/o-^2c<;Mod / 

Indicating with -Ri^Modia^Mag noise spectral density at the signal frequencies, and as- 
suming i2<^Mod+2f^Mag = ^a;Mod-2nMag) ^^c seusitivlty spectral density ^'sens of the ellipsometer 
for a unity signal to noise ratio is 

D 

,T, -"'<^Mod+2nMag /'00\ 

*Sens = , {28} 
V^-'0-'2a;Mod 

In principle the r.m.s. noise limit for such a system is determined by the r.m.s. shot- 
noise ishot of the DC current 2dc generated by the modulation amplitude by the 
extinction ratio loqcr"^ and the by DC component of the spurious ellipticity I^qa^Q (see table 

I 1 

Zshot = V2e2DcAz/ = y 2e/og(a2 + ^ + «?)c) (29) 

where g is the quantum efficiency of the photodetector, Ai/ is the bandwidth and e is the 
electron charge. In the case Q ^ cr^ and Q ^ o^q the DC current will only depend on C,q 



and by substituting -R^Modif^Mag ~ "^shot/ (^v Ai/) into equation (28) the shot-noise sensitivity 
spectral density \E'shot becomes 



For a power Jq = 10 mW and a quantum efficiency g = 0.7 A/W this leads to a sensitivity 
spectral density of \E'shot — 3.4 • 10^^;^^. It is interesting to note that such a limit depends 
exclusively on the laser power before the analyzer and the quantum efficiency of the detector. 

Other intrinsic noise sources are photodiode dark current noise idark = VdarkV Az//G, 
Johnson current noise = ^JAKbTAv /G of the transimpedence G in the amplifier of the 
photodiode, and residual laser intensity current noise zrin = IqQ ■ RlN{uj)\^ Av. These noises 
must be kept below ishot at a frequency near c^Mod in order to reach the theoretical sensitivity. 
The expressions for these noise contributions to the ellipticity spectral noise density can be 
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obtained from equation (28): 



*shot 
^dark 




2KbT 1 



RIN 



G /ogCo 
RIN(u;Mod) V'(^' + CoV2)' + (CoV2)' 



(31) 
(32) 
(33) 

(34) 



72 Co 

With the PVLAS experimental parameters given in Table [IT], the contribution of each of 
these noises to the sensitivity spectral density can be plotted as a function of the modulation 
amplitude Co- This allows the optimization of the modulation amplitude. Figure [5] shows 
the corresponding plots for the infrared and green configurations. In each graph a cross 
marks the current experimental sensitivity. 

Table II: Experimental parameters for the two laser configurations. 



Configuration 


Photodiode 


Cavity output 


Extinction 




Gain G 


Photodiode noise 




efficiency q [A/W] 


power Iq [mW] 


ratio cr^ 


[1/VHz] 


[V/A] 




Green 


0.2 


1.5 


5 • 10-7 


2 • 10-5 


109 


2 


infrared 


0.7 


60 


5 • 10-7 


2 • 10-5 


107 


8 



As can be seen, in both configurations we are still well away from the theoretical limit. 
The one noise source which cannot, at the moment, be controlled is the low frequency 



spurious ellipticity a{t) (see equation (22)) induced in the system. We believe this noise is 
due to the movement of the granite tower. Since the mirrors and optical elements have a 
structural birefringence 'map' with a gradient [20], these movements will generate ellipticity 
noise. Indeed, we have measured the induced ellipticity as a function of movement at the 
top of the tower and found a value ~ 0.4 m~^. To reach the theoretical limit, the relative 
movement of the top of the tower respect to the lower optical bench must be less than 
10-« m/VRi. 
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Figure 5: Calculated noise contributions for the infrared and green configurations of the PVLAS 
apparatus. See text and table\l^for the parameters of each configuration. 



III. PREVIOUS PHOTON-PHOTON SCATTERING RESULTS 

In this section we review the resuhs obtained by D. Bernard et al. in a direct search for 
diffused photons from two colhding beams [9] and the results from the BFRT collaboration 
[21], another experiment designed to search for axions via ellipsometric techniques. There 
are several proposals for improving the direct scattering measurements [221 l23l [211 [25] and 
for detecting the magnetic birefringence of vacuum [251 1221 12H] • 



A. Stimulated photon-photon scattering 

In [0] direct photon-photon scattering was searched for. Differently from what was pre- 
viously done [TD], the authors searched for stimulated scattering when three high power 
ultra short beams were crossed. The advantage of this configuration is to fix the angle 
and wavelength for the scattered photon: scattered photons are searched for in a defined 
direction and at a defined wavelength. Indeed in two-beam scattering, energy (indicated 
as Cj) and 3-momentum (indicated as kj) must be conserved: ei + 62 = 63 + 64; ki + k2 
= ks + k4. This leaves two free parameters for the final state. In the three beam con- 
figuration the fourth (scattered) beam must satisfy the condition k4=ki + k2 - ka and 
A4 = (1/Ai + I/A2 — I/A3) \ In the configuration being discussed Ai = A2 = 800 nm and 
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As = 1300 nm, resulting in A4 = 577 nm. 

The theoretical analysis of three wave mixing in vacuum develops in analogy with three 
wave mixing in a medium [H]. In this latter case the medium polarizability is written as 

V{t) = x}'^£it) + f^£\t) + f^S^it) + ... (35) 

In four wave mixing the authors show that the growth rate of the electric field Sq4 of the 
scattered beam depends on x^^^ and can be written as 

dSo4 (3)c c c -.1, d did 

When considering the Euler-Heisenberg Lagrangian density correction one finds that the 



growth rate of £^04 in vacuum has the same form as equation (36 ) with the QED susceptibility 
of vacuum having the expression 

X^'^ = XJ'^ = .JtTr = - 3.0 ■ lO-^'K mVV^ (37) 

where is a parameter depending on the direction of the incident beams and their polar- 
ization. In the experiment reported in [22] K ^ 0.56. 



Integration of equation (36) leads to a number of counts per pulse crossing 



= ^PMes,e,oss—,^,^^^,^{x^ ^) (38) 

where Ei are the energies of the three incoming laser pulses, epM,^Sp, ^loss are the quantum 
efficiency of the photomultiplier tube, the transmission of the spectrometer and a loss factor 
due to a beam position oscillation, and w and cr are respectively the beam waist and 
bunch length. The value of the third order susceptibility x^^^ of nitrogen was measured and 
compared to other experiments. An order of magnitude agreement was observed, allowing a 
calibration of the apparatus. A comparison between the expected vacuum counts calculated 
from QED and the observed counts resulted in a limit on the total photon-photon cross 
section at 0.8 eV center of mass energy of 



mernard) ^ ^Hob^ fQED) ^ ^ 5 . ^2 (39) 

A^4,QED 



which is eighteen orders of magnitude larger than the theoretical QED cross section. 
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B. Brookhaven-Fermilab- Rochester- Trieste (BFRT) results 



The principle of the ellipsometer in the BFRT collaboration is the same as the one shown 
in Figure |2j In this case the cavity was a multi-pass cavity with a number of reflections which 
varied from 34 to 578. The laser wavelength was 514.5 nm and the length of the magnetic 
field region was 8.8 m. To modulate the magnetic vacuum birefringence the magnetic field 
was ramped from 2.63 to 3.87 T at a frequency of 30 mHz. The sensitivity of the apparatus 
varied as a function of the number of refiections in the multi-pass cavity and consequently 



did the final limit on the acquired ellipticity. The results are summarized in table III 



Table III: Summary of the BFRT experimental parameters and results together with the limit 
achieved on the parameter Ag. 



Number of 


Measured 


Ellipticity upper 


Ae Upper 


passes 


sensitivity [l/-v/Hz] 


bound V'limit at 95%C.L. 


bound [r-2] 


(shunt) 


2.6 • 10-*^ 


7.7-10-10 


n.a. 


34 


7.9 • 10-*^ 


2.0-10-9 


1.4 - 10-19 


578 


1.5 • 10-6 


5.1 - 10-8 


2.1 - 10-19 



Due to the fact that the magnetic field was ramped around a central value Bq = 3.25 T, 
with an excursion ±AB = ±0.62 T, the expression for the Euler-Heisenberg induced mag- 
netic birefringence is 

An = 3Ag{2BoAB) (40) 



and the limit on A^. attainable from the BFRT results is 

A 



limit ' 



(41) 



'GBqABNL 

When translated into a photon-photon cross section for unpolarized light these results 
give a limit of 

a™ < 1.6 ■ 10-^^ m^ (42) 

This must be compared to the QED photon-photon cross section at the same wavelength of 
514.5 nm which is 4?™^ = 1.44 ■ 10-^^ m^. 
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IV. PVLAS RESULTS 



Gas measurements, for calibration, and vacuum birefringence measurements were con- 
ducted with the apparatus in both the infrared and green configurations. We present here 
measurements taken with the magnet energized at 2.3 T. This choice of field strength is 
motivated by the strong suppression of the stray field outside the magnet. Indeed, at higher 
fields the presence of a stray field has resulted in a yet to be understood spurious ellipticity 
signal [15] . The total integration times were Tio64 = 45200 s at 1064 nm and T532 = 28300 s 
at 532 nm. 

A. Gas calibration measurements 

Calibration of the ellipsometer is done by taking advantage of the Cotton-Mouton effect 
[16] in gases. In the presence of an external magnetic field perpendicular to the propagation 
of a light beam, gases become birefringent. Depending on the gas, the induced birefringence 
may be positive (ny — n_L > 0; e.g. He) or negative {n^\ — n± < 0; e.g. N2). These 
measurements also allow the verification of the Fourier phase of the sidebands of c^Mod at 
c^Mod ± Sf^Mag with what they should be (see table [l]). Indeed, the ellipticity induced by a 
birefringence is maximum when the angle between the polarization and the slow axis defined 
by the magnetic field is 45° (n|| — > 0). Figure [6] shows a polar plot corresponding to the 
amplitude and phase of the signal demodulated at c^Mod due to Helium gas at four different 
pressures (5, 10, 15 and 20 mbar), measured with a field intensity of 2.3 T. As can be seen, 
the experimental values lie on a straight line with a Fourier phase of 125°. This defines the 
physical axis of signals. A gas with a negative Cotton-Mouton constant would generate a 
signal at 180° with respect to the signals shown in Figure [6j 

Different gases were measured [T71 [18] resulting in an accuracy better than 20%. 

B. Vacuum measurements 

The complete data sets of the signals from the lock-in amplifiers demodulated at c^Mod 
were analysed by a Fourier transform. No peak was found at 2Quag as would be expected 
from a magnetically induced ellipticity. The data will therefore be presented as a noise his- 
togram in a frequency band around 2nMag, between 1.92 l^Mag and 2.08 ^Mag (see Figure [?]). 
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Figure 6: Polar plot for the ellipticity signal generated with a 2.3 T magnetic field intensity when 
Helium gas is present in the vacuum chamber. The figure shows the signal for four different gas 
pressures: 5, 10, 15 and 20 mbar. Each data point represents amplitude and phase of the signal 
peak observed in a 100 s long time record. 

The probability density function for the Fourier amphtude rp = \/xp + t/p, where xp and 
y-p are the projections of the Fourier transform along the physical and the non physical phases 
(defined above) respectively, is given by the Rayleigh distribution, pi^rp) = r-pe~^ /a'^. In 
this expression a is the standard deviation of the Gaussian noise distributions of x-p and 
Uy from which we deduced our limits on the induced ellipticity. To extract a from the 
histograms the r-p noise distributions were fitted with the Rayleigh distribution. 

The 95% confidence limits are then deduced from the cumulative distribution P(rF) = 
1 — e~''F/^°'^ . In Figure [t] are shown the histograms and fits for the measurements taken with 
the 1064 nm and 532 nm lasers. Superimposed on these graphs, represented by a vertical 
black line, are the values at the bins corresponding to 2VLMag of the Fourier spectrum of 
the demodulated signal. As can be seen, these are well within the noise distributions. The 
standard deviations of the two distributions are very similar even though the integration 
time with the 532 nm laser was 70% of the integration time with the 1064 nm. This is due 
to the lower noise encountered with the 532 nm setup. 



Table IV gives the 95% confidence level background values for the ellipticity measurements 



with the PVLAS and the BFRT apparatus. The stimulated scattering results are also 
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Ellipticity Ellipticity 

Figure 7: Noise distributions in the magnet rotation frequency hand 1.92 ^Mag - 2.08 ^Mag for the 
2.3 T ellipticity measurements with the IR (left) laser and the green (right) laser. The vertical line 
indicates the value in the Fourier spectrum corresponding to 20, Mag- Indicated in the legend is the 
value of a for the two wavelengths. 



Measurement 
type 


Photon 
energy 


noise 
floor 


(i? ^)cquiv 

[T^m] 


Ae bounds 
95% C.L. [T-2] 


0"^^ bounds 
[m2] 


/ (QED) 


stimulated scatter [9] 
BFRT Ellipticity [H] 
PVLAS Ellipticity [E] 
PVLAS Ellipticity 


0.8 eV c.o.m 
2.42 eV 
1.17 eV 
2.34 eV 


4.9 • 10-9 
1.4 • 10-^ 
1.4 • 10-s 


1197 
238000 
124000 


1.2 • 10-15 
1.4 • 10-19 
6.6 • 10-21 

6.3 • 10-21 


1.5 • 10-52 

1.6 • 10-57 

4.6 • 10-62 

2.7 • 10-60 


8 • lOl'^ 
11 • 109 
2.5 • 10^ 
2.3- 10^ 



Table IV: Ellipticity results and stimulated scatter results for a. 



reported where appropriate. The parameters for the different configurations are also reported 
in the same table. In the last three columns we also report the limits on Ag, a^^ and the 
ratio o^^l o^f''^^ . 

Although experimental results have not reached the predicted QED values, bounds have 
been improved. We believe that at the moment the sensitivity is limited by seismically 
induced spurious ellipticities. 



V. DISCUSSION AND CONCLUSIONS 



We have reported here the interpretation of vacuum magnetic birefringence limits in terms 
of photon-photon scattering. Although the sensitivity of our apparatus has not reached its 
theoretical shot noise limit, the ellipso metric technique is at the moment the most sensitive 
one for approaching low energy non linear QED effects. In a general post-Maxwellian frame- 
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work, though, direct scattering measurements are necessary to extract the free parameters 



^rji and ^772 (see equation ( 10 )). In the Euler-Heisenberg framework we are now a factor 4800 
away from the theoretical parameter, = 1.32 • 10^^^ T~^, describing non hnear quantum 
electrodynamic effects: 

< 6.3 • 10"^^ T-^ @ 95% C.L. (43) 

Always in the Euler-Heisenberg framework, from the experimental bound on Ae one can 
place the following upper bounds on the photon-photon cross section for non polarized light 
in the limit hu <^ rrieC^, at 1064 nm and 532 nm respectively of: 

CT^!;^^^) < 4.6 ■ 10-62 ^2 (44) 
^(532) < 2.7 ■ 10-*^° m^ (45) 
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